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SECTION  1 . INTRODUCTION 


In  a recent  paper  (Greifinger  and  Greifinger  [1978] , 
hereinafter  referred  to  as  I)  , the  authors  derived  an  approx- 
imate expression  for  the  complex  propagation  constant  for  ELF 
propagation  in  the  Earth-ionosphere  waveguide.  The  theory 
was  developed  for  an  ionosphere  that  was  sufficiently 
disturbed  that  there  was  no  significant  penetration  of  the 
electromagnetic  field  to  altitudes  where  anisotropy  due  to 
the  Earth's  magnetic  field  had  to  be  taken  into  account.  It 
was  shown  that  the  propagation  constant  is  determined  by  four 
parameters,  two  frequency-dependent  altitudes  and  the  local 
conductivity  scale  heights  at  those  altitudes.  The  lower 
altitude,  denoted  by  hQ,  is  where  the  conduction  current 
becomes  equal  to  the  displacement  current,  and  the  higher 
altitude,  denoted  by  h^,  is  where  the  local  reciprocal  wave 
number  becomes  equal  to  the  local  scale  height  of  the  refrac- 
tive index.  There  is  a Joule  heating  maximum  at  hQ  arising 
from  predominantly  vertical  currents,  and  a secondary  maximum 
in  the  vicinity  of  h^  arising  from  predominantly  horizontal 
currents.  At  altitudes  in  the  vicinity  of  and  below  hQ,  the 
electric  field  is  predominantly  vertical,  and  energy  flow  in 
the  waveguide  is  in  the  horizontal  direction.  Within  a few 
scale  heights  above  hQ,  the  electric  field  becomes  horizon- 
tally polarized,  and  energy  flow  is  vertical. 

In  this  paper,  the  theory  developed  in  I is  extended  to  include 
the  effects  of  anisotropy.  For  the  sake  of  mathematical 
simplicity,  it  is  assumed  that  the  Earth's  field  is  vertical, 
which  limits  the  validity  of  the  theory,  strictly  speaking, 
to  high  magnetic  latitudes.  The  generalization  to  arbitrary 
dip  angle  is  straightforward,  though  somewhat  tedious.  It 
turns  out,  however,  that  there  is  no  significant  dependence 
of  the  propagation  constant  on  magnetic  latitude  except  very 
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close  to  the  magnetic  equator,  so  the  results  have  a rather 
large  geographic  range  of  applicability. 

When  anisotropy  is  included,  the  parameters  which  determine 
the  propagation  constant  differ  for  daytime  and  nighttime 
ionospheric  conditions.  In  both  cases,  two  parameters  which 
enter  are  the  frequency-dependent  altitude  hQ  at  which  the 
conduction  current  parallel  to  the  magnetic  field  becomes 
equal  to  the  displacement  current,  and  the  local  scale  height 
of  the  parallel  conductivity  a . Under  daytime  conditions, 
two  additional  pairs  of  altitudes  (frequency-dependent)  and 
scale  heights  appear  as  parameters.  One  pair  is  the  altitude 
at  which  the  local  reciprocal  wave  number  for  vertically 
propagating  0 waves  becomes  equal  to  the  local  scale  height  of 
the  refractive  index,  and  the  scale  height  of  this  refractive 
index.  The  other  pair  are  the  corresponding  quantities  for 
vertically  propagating  X (whistler)  waves.  It  is  assumed  in 
this  paper  that  these  altitudes  are  attained  in  a region  of 
the  ionosphere  where  l°Hl  >5>  I ap  I ' aH  bein9  the  Hall  conduc- 
tivity and  od  the  Pedersen  conductivity.  Under  these  condi- 
tions, which  apply  over  a substantial  altitude  range,  the  two 
pairs  of  parameters  become  identical  and  the  analysis  is  some- 
what simpler.  The  vicinity  of  the  single  altitude  h^  is  in 
this  case  a region  of  reflection,  rather  than  of  significant 
heating,  as  was  the  case  for  the  isotropic  ionosphere.  The 
0 wave  undergoes  nearly  total  reflection  at  this  altitude, 
but  some  of  the  X wave  energy  leaks  out  of  the  waveguide, 
thereby  contributing  to  the  attenuation. 

For  typical  nighttime  conditions,  a sharp  reflecting  E-region 
bottom  may  be  encountered  before  the  local  reciprocal  wave 
number  becomes  equal  to  the  local  scale  height.  Under  such 
circumstances,  the  altitude  of  the  E-region  bottom  replaces 


as  a parameter,  and  the  refractive  index  scale  height  is 
replaced  by  the  local  wavelength  on  the  E-region  side  of  the 
bottom. 

As  in  I,  the  approximate  propagation  constants  are  given  by 
simple  algebraic  expressions  involving  the  various  parameters. 
Numerical  results  have  been  obtained  for  a hypothetical  day- 
time ionospheric  profile  in  which  the  electron  density  is 
assumed  to  increase  exponentially  with  altitude  and  the  elec- 
tron collision  frequency  to  decrease  exponentially  with 
altitude.  The  results  are  in  excellent  agreement  with  full- 
wave  calculations  for  the  same  profile  which  were  carried  out 
by  Dr.  William  Moler  of  Naval  Ocean  Systems  Center.  The 
agreement  lends  support  to  the  physical  assumptions  on  which 
the  theory  is  based. 

Booker  and  Lefeuvre  (1977)  have  also  proposed  a method  for  the 
calculation  of  approximate  ELF  propagation  constants  in  the 
anisotropic  Earth-ionosphere  waveguide.  The  method  is  based 
on  a simplified  model  of  the  waveguide  in  which  the  ionosphere 
is  cut  off  discontinuously  at  a frequency-dependent  level  and 
abolished  below  this  level.  The  altitude  at  which  the  ion- 
osphere is  truncated  corresponds  very  closely  to  the  altitude 
h-^,  and  their  treatment  of  the  fields  above  this  level  is 
quite  similar  to  ours.  In  their  formulation,  the  region 
around  hQ,  which  plays  a unique  and  important  role  in  our 
theory,  is  treated  as  part  of  the  vacuum.  This  results  in  a 
predicted  phase  velocity  which  is  in  general  significantly 
higher  than  that  obtained  from  our  theory.  They  attempt  to 
account  for  the  attenuation  due  to  Joule  heating  at  levels 
below  their  truncated  ionosphere,  but  do  so  in  a manner  which 
may  seriously  overestimate  the  size  of  the  effect,  especially 
at  the  low  end  of  the  ELF  band. 
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SECTION  2.  BASIC  EQUATIONS 

It  will  be  assumed  that  the  Earth-ionosphere  waveguide  is 
horizontally  stratified  and,  as  discussed  in  the  Introduction, 
that  the  geomagnetic  gield  is  in  the  vertical  (z)  direction. 

The  equations  governing  electromagnetic  propagation  in  the 
anisotropic  waveguide  are,  of  course,  Maxwell's  equations 

^ X g = (1) 

v x 5 = - i-*$-  e , (2) 

c 

where  we  have  assumed  a time  dependence  e , and  the 
generalized  Ohm's  law 

? = a0&  * + ap&  “ * W + aH(^  x ez>  ' (3) 

where  §_  is  a unit  vector  in  the  z direction  and  o , o , and 
z op 

a„  are  the  parallel,  Pedersen,  and  Hall  conductivities, 
respectively.  The  fields  can  be  written  in  terms  of  the  cus- 
tomary potentials  as 


1 

S = ? x A 

(4) 

2 = ^ + io»X  . 

(5) 

Combining  Eq.  (4)  with 

the  time  derivative  of  Eq. 

(2) , we 

obtain 

i(j[V(^»A)  - 

V2A]  - iuQa»l  ~ k2^  = 0 . 

(6) 

I 

! 


The  modal  solutions  of  this  equation  are  obtained  by  a 
separation  of  variables.  Th-2  rectangular  components  of  A and 
the  potential  ^ satisfy  a two-dimensional  wave  equation 

(V2  + k2S2)F  = 0 , (7) 

2 

where  is  the  Laplacian  operator  in  the  horizontal  plane  and 
kS  is  the  (complex)  horizontal  propagation  constant  for  the 
mode . 


At  ELF,  only  the  lowest  mode  is  non-evanescent  in  the  Earth- 
ionosphere  waveguide.  In  the  absence  of  anisotropy,  this  is  a 
TM  mode  which  is  derivable  from  a vector  potential  with  only 
a vertical  component  A . When  anisotropy  is  included,  the 
lowest  mode  acquires  a TE  component,  which  requires  a horizon- 
tal component  of  the  vector  potential.  It  will  be  shown  that 
all  the  boundary  conditions  can  be  satisfied  by  a vector 
potential  of  the  form 


A = Ae  - 
z z 


TW  * x {UV 


where  u is  a scalar  function. 


(8) 


The  various  relations  above  may  now  be  incorporated  into 
Eq.  (6)  to  obtain  a set  of  coupled  equations  for  the  system. 
The  vertical  component  of  Eq.  (6)  becomes 


3+  s. 
3z 


(9) 


while  the  horizontal  component  becomes 
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(S  x ^)  [V2u  + (iu  wo  +k2)u  - iu  u>au<J)] 

Z Op  Orl 

- $L  [iw$  • A - (iyowop  + k2)^  - iyocooHu]  = 0 . 

(10) 

A gauge  condition  remains  to  be  specified.  A convenient  choice 
is  clearly 


iu>^  • A = (ip  wo  + k )il)  + iy  ii)o..u  , (11) 

op  o H 

which  removes  the  terms  involving  from  Eq.  (10).  The 
contents  of  the  first  brackets  must  then  also  vanish,  which 
gives 


— 7 = ~ [iy_.too_  + k2  (1  - S2)  ] u + iu  o>o„i|)  . (12) 

3z2  p ° H 

The  basic  equations  are  (9)  , (11) , and  (12)  , in  which  S 
appears  as  an  eigenvalue  parameter.  The  eigenvalue  is  deter- 
mined by  solving  these  coupled  equations  with  appropriate 
boundary  conditions  at  the  ground  and  at  large  altitudes.  An 
approximate  method  for  determining  the  eigenvalue,  similar  to 
that  developed  in  I,  will  be  presented  below.  The  method 
consists  of  constructing  two  approximate  analytic  solutions 
of  the  equations.  One  solution  obeys  the  proper  boundary 
conditions  at  the  ground  and  is  valid  up  to  an  altitude  a few 
scale  heights  above  hQ.  The  other  obeys  the  proper  boundary 
conditions  at  large  altitudes  and  is  valid  down  to  an 
altitude  a few  scale  heights  below  h^.  There  is  a common 
altitude  region  where  both  solutions  are  valid.  The  eigen- 
value is  determined  by  requiring  that  the  leading  terms  of  the 
two  solutions  agree  in  the  overlap  region. 
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SECTION  3.  APPROXIMATE  SOLUTIONS  AND  EIGENVALUES 

3.1  VACUUM  AND  LOWER  IONOSPHERE 


The  conductivities  of  the  medium  are  increasing  functions  of 
altitude.  As  in  I,  we  define  a frequency-dependent  altitude 
hQ  at  which  the  parallel  conduction  current  becomes  equal  to 
the  displacement  current.  Between  the  ground  and  a few 
conductivity  scale  heights  above  hQ,  the  right-hand  sides  of 
Eq.  (11)  and  (12)  are  very  small,  and  these  equations  become 
approximately 


(13) 


- 2 
3 u 

: 2 

dZ 


= 0 


(14) 


At  ELF,  the  ground  may  be  considered  as  perfectly  conducting, 

and  the  appropriate  boundary  condition  is  therefore  E ,.(0)  =0. 

This  in  turn  requires  iMO)  = u(0)  = 0.  The  solution  of 

Eqs . (9)  , (13)  , and  (14)  satisfying  these  boundary  conditions 

is 


A = 
z 

A (constant) 

o 

(15) 

u = 

Bz 

(16) 

ip  = 

-iuiA 

o 

z - s 2 fZ 

dz 

(17) 

° 1 

r co“/ 
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The  integrand  in  Eg.  (17)  is  essentially  unity  up  to  a 
conductivity  scale  height  or  so  below  h^  and  becomes  very 
small  a few  scale  heights  above  this  altitude.  It  is  there- 
fore necessary  to  represent  oQ  accurately  only  in  the  neigh- 
borhood of  hQ,  where  a simple  exponential  provides  a very  good 
approximation.  Thus,  we  write 


a = e toe 
o o 


(z'ho)Ao 


where  is  the  local  scale  height  at  hQ.  The  integral  in 
Eg.  (17)  can  now  be  evaluated,  giving 


-iwA  [z 
°( 


S2  h 

o 2 ^o 


/ 

S0ln  ^1-ie 


(^-2)  A, 


The  constant  £ and  the  eigenvalue  S remain  to  be  determined. 


3.2  MIDDLE  AND  UPPER  IONOSPHERE 


A few  scale  heights  above  hQ,  the  parallel  current  becomes 
much  larger  than  the  displacement  current,  and  Eg.  (9) 
becomes  approximately 


3 

a7  + 1U)AZ  = 0 <2°: 

The  left-hand  side  of  Eg.  (20)  is  exactly  the  vertical  compo- 
nent of  the  electric  field,  which  thus  becomes  very  small 
somewhat  above  hQ.  The  electric  field  thus  undergoes  an 
important  transition  from  vertical  to  horizontal  polarization 
in  the  region  around  hQ. 


At  the  altitudes  in  question,  the  Pedersen  and  Hall  currents 
are  also  much  larger  than  the  displacement  current,  so  that 
the  terms  proportional  to  k in  Eqs . (11)  and  (12)  may  be 
neglected  (a  QL  approximation) . These  equations  therefore 
become 


yf  + iwo“V  + itJoWaHU  = 0 


~2  + iyoWapU  " = 0 ' 

o z 


where  A has  been  eliminated  by  use  of  Eq.  (20). 


If  we  introduce  as  variables 


ip+  = ip  ± iu  , 


we  obtain  the  two  uncoupled  wave  equations 


a2^+ 


+ 2 2 
2~  + n+k  Y + = 0 


where 


n±  = ‘ "h' 

O 


The  upper  sign  corresponds  to  vertical  0 wave  propagation  and 
the  lower  sign  to  vertical  X wave  (whistler)  propagation . 
These  equations  must  be  solved  subject  to  the  radiation 
boundary  condition  at  large  altitudes. 
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3.2.1  Daytime  Ionospheric  Conditions 


Following  the  procedure  of  I,  we  will  construct  an  approximate 
solution  of  the  upper  ionosphere  equations  which  satisfies  the 
boundary  condition  at  large  altitudes  and  which  has  a common 
region  of  validity  with  the  lower  ionosphere  solution.  This 


involves  the  introduction  of  the  altitude  h^,  at  which  the 


local  wave  number  |n|k  is  equal  to  the  reciprocal  of  the 
refractive  index  scale  height.  Under  highly  disturbed  daytime 
conditions,  this  occurs  at  an  altitude  where  the  ionosphere  is 
nearly  isotropic,  i.e.,  where 


2 2 io 

n,  s n = 

+ - e cj 

o 


This  was  the  case  treated  in  I,  for  which  the  eigenvalue  was 


shown  to  be 


s2  . (hl  * ^l) 


(ho  * ❖ „) 


where  is  the  conductivity  scale  height  at  the  altitude  h^. 


The  altitudes  h and  h.  were  shown  to  be  locations  of  maximum 
o 1 


Joule  heating,  the  lower  altitude  maximum  being  associated  with 
vertical  currents  and  the  higher  altitude  maximum  with  horizon- 
tal currents.  There  was  no  significant  field  penetration  much 


above  the  altitude  h^. 


Under  normal  daytime  conditions,  and  over  most  of  the  ELF  band. 


the  altitude  h-^  occurs  where  >>  |Op|,  i.e.,  where 


2 2 _ h ; 

n s -n  = = -n 

+ - ew 

o 


I 

( I 


1 

In  the  vicinity  of  h^,  we  may  approximate  the  index  of 
refraction  by  an  exponential  with  a scale  height  appropriate 
to  that  altitude.  Thus  we  write 

2 2 <z-hi)/^ 

n = n^e 

(29) 

where  is  defined  as  the  frequency-dependent  altitude  at 

which 

2|n1|k?1  = 1 . 

(30) 

i 

(The  scale  height  for  the  refractive  index  has  been  taken 
2^1  for  consistency  with  the  theory  for  the  isotropic  ion- 
osphere . ) 

as 

i 

Since  the  imaginary  part  of  oH  is  very  small,  the  quantity  n.^  is 
nearly  real.  The  outgoing  wave  solutions  of  Eq.  (24)  are  then 

ip+  = aH(J1)  (iy) 

(31) 

i 

ip_  = yaH^1)  (y) 

(32) 

i 

i 

(z-h.)/2;1 

y = e 

(33) 

where  a and  y are  constants  to  be  determined.  The  potential 
functions  ip  and  u,  as  determined  from  Eq.  (23),  are 


* ’ §[' 

(iy)  + YHq1)  (y)] 

(34) 

u - A 

[«oX)  (iy)  - (y)j  . 

(35) 
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The  potentials  given  by  Eqs.  (34)  and  (35)  must  agree  with 
Eqs.  (16)  and  (19),  respectively,  in  the  altitude  range 
hQ  <<  z <<  h^,  where  both  approximations  are  valid.  In  this 
altitude  range,  |y|  <<  1 and  the  Hankel  functions  may  be 
approximated  by  their  small  argument  expansions.  The  leading 
terms  give 


u — [(1  + y)  (z  - hi>  + i-TTYsJ 

*-2%[(1-Y)(z-hl)  ' iwl] 

(hQ  <<  z <<  h-j^)  . 


(36a) 


(36b) 


Matching  these  functions  to  the  leading  terms  of  the  lower 
ionosphere  solutions  in  the  same  altitude  range,  we  obtain 


Y 


+ nrcx) 


(37) 


S2  = 


hl(hl  + ilT?l) 


(ho  “ (hl  + ¥*!> 


(38) 


Although  this  completes  the  derivation  of  the  approximate 
eigenvalue,  there  are  two  additional  relationships  which  exist 
among  the  three  remaining  undetermined  constants.  These  are 


a 


C.(h.  + iTre.) 

2"“  Rh1  +~rnjr 


(39) 
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(2h1  + iir 


with  A remaining  as  an  arbitrary  normalization  factor.  An 
approximate  analytic  solution  has  thus  been  obtained  not  only 
for  the  eigenvalue,  but  for  the  height  dependence  of  the  var- 
ious field  components  as  well. 

The  approximate  eigenvalue  given  by  Eq.  (38)  is  slightly  dif- 
ferent from  its  counterpart  for  the  isotropic  ionosphere  given 
by  Eq.  (27).  However,  since  C0/hQ  <<  1 and  <<  both 

eigenvalues  are  to  first  approximation 


S 


iTT  l^O 


(41) 


The  relative  phase  velocity  is  thus  approximately 


(42) 


and  the  horizontal  attenuation  rate  in  decibels  per  megameter 
is  approximately 


a 


(43) 


where  f is  the  frequency  in  Hertz.  Thus,  in  first  approxima- 
tion, the  relative  phase  velocity  depends  only  on  the  ratio 
of  the  two  altitudes,  and  is  independent  of  the  scale  heights. 
This  is  quite  different  from  the  theory  of  Booker  and  Lefeuvre 


(1977),  in  which  the  relative  phase  velocity  differs  from  unity 
by  an  amount  proportional  to  C-^/h^. 

The  horizontal  attenuation  in  the  waveguide  arises  from  two 
sources.  The  part  proportional  to  C0/hQ  is  due  to  Joule 
heating  by  vertical  currents,  which  has  a maximum  at  hQ  as  in 
the  isotropic  case.  However,  the  part  proportional  to 
is  no  longer  associated  with  Joule  heating  by  horizontal  cur- 
rents. The  horizontal  currents  are  now  Hall  currents,  which 
are  non-dissipative . The  field  at  this  altitude  has  been 
decomposed  by  the  anisotropic  medium  into  two  vertically  prop- 
agating waves,  an  0 wave  component  whose  index  of  refraction 
is  almost  purely  imaginary  and  an  X wave  component  whose  index 
is  almost  purely  real.  Thus,  the  0 wave  undergoes  total 
reflection  in  this  region,  whereas  some  of  the  X wave  energy 
leaks  out  of  the  waveguide.  The  part  of  the  attenuation  rate 
proportional  to  ^/h^  is  associated  with  this  leakage,  as 
pointed  out  by  Booker  and  Lefeuvre  (1977) , who  obtained  a very 
similar  result  for  this  quantity.  The  reflected  part  of  the 
X wave  results  in  a small  TE  component  in  the  field  at  the 
ground.  From  Eqs.  (4) , (8) , (16) , and  (40) , the  ratio  of  the 
TE  and  TM  components  at  the  ground  is 

^TE  = / 1_  9u\  B _ , ^1  (44) 

BTM  l ia)Az  8z/z=0  ia)Ao  (2hi  + i^i) 

The  magnetic  field  at  the  ground  thus  has  slight  elliptical 
polarization,  with  the  major  axis  nearly  perpendicular  to  the 
plane  of  propagation. 

The  theory  has  been  developed  under  the  assumption  that  h-^  is 
reached  at  an  altitude  where  |o_J>>|a  I.  If  the  criterion  is 
met  where  the  Pedersen  conductivity  is  not  negligible,  there 


is  a separate  pair  of  values  of  and  for  the  0 and  the  X 
components,  which  in  general  differ  only  slightly.  It  is  not 
difficult  to  generalize  the  theory  to  this  situation.  This 
results,  however,  in  only  a small  correction  to  the  eigenvalue 
obtained  by  replacing  the  two  pairs  of  parameters  by  their 
average. 

3.2.2  Nighttime  Ionospheric  Conditions 

Under  nighttime  conditions,  a sharp  E-region  bottom  is  usually 
encountered  before  the  altitude  h-^  is  established.  The  elec- 
tron density  undergoes  a very  sharp  increase  in  passing  through 
the  bottom,  above  which  it  can  be  quite  variable.  We  will 
consider  only  the  simplest  model  where  the  density  above  this 
bottom  varies  slowly  on  the  scale  of  the  local  wavelength. 

Under  such  conditions,  the  phase  integral  approximation  is 
valid  above  the  E-region  bottom,  and  we  may  take  as  solutions 
of  Eq.  (24) 


= 


= A 


1 ( r z 

+ (n  k)  2 exp'ik  / n (z)dz 

''  ’ ' I JhE  ‘ 


(45) 


where  h„  is  the  altitude  of  the  E-region  bottom  and  A are 
£ 

constants.  The  associated  potential  functions  <p  and  u have 
the  form 


\p  = A(^_  + 6^+) 


u = iA(^_  - 6i(>+) 


(46) 


(47) 


where  A and  6 are  constants  to  be  determined.  Under  nighttime 
conditions,  the  solution  up  to  the  E-region  bottom  is  well 
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approximated  by  Eqs.  (16)  and  (19) . The  boundary  conditions 
require  that  the  tangential  electric  and  magnetic  fields  be 
continuous  at  z = h 


E ' 


. du 
of  u'  dP 


This,  in  turn,  requires  the  continuity 

If  we  assume,  as  in 


and  ^jr  at  this  interface. 


Eq.  (28),  that 


n+(hE)  S n^(hE)  * 


-nr 


(48) 


then  application  of  the  boundary  conditions  leads  to 


S 


2 


+ ! 


1 + i(l+2e)  |) 
1 + |e(l+i)|) 


where 


kVE  * 


(49) 


(50) 


Assuming  that  n_  is  nearly  real,  the  eigenvalue  is  in  first 
approximation 


1 

4knEhE 


)] 


(51) 


Comparison  with  the  daytime  results  shows  that  the  altitude  of 
the  E-region  bottom  has  replaced  the  frequency-dependent 
altitude  h^  as  a parameter  and  the  local  wavelength  just  inside 
the  E-region  has  replaced  The  physical  processes  occurring 

at  hE  are  the  same  as  those  associated  with  h^,  namely  nearly 
total  reflection  of  the  0 wave  and  partial  leakage  of  the 
X wave.  Although  we  have  considered  only  the  simplest  case  of  a 
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SECTION  4.  NUMERICAL  RESULTS  AND  CONCLUSIONS 


The  validity  of  the  theory  has  been  examined  by  comparing 
approximate  eigenvalues  with  a full-wave  calculation  for  a 
hypothetical  daytime  conductivity  profile.  For  simplicity, 
it  was  assumed  that  only  electrons  contributed  to  the 
conductivity  at  all  altitudes.  Although  this  assumption  may 
be  unrealistic,  it  does  not  invalidate  a comparison  of  an 
approximate  solution  with  a full-wave  solution  for  the  same 
profile.  The  electron  density  and  collision  frequency  pro- 
files were  assumed  to  be  simple  exponentials 


zAe 

Noe  ' 

(52) 

-zA 

V e v 

o 

(53) 

with  N = 3.73  x 10^/m"^,  v = 1.63  x lO^sec  \ and  £ = £ 

o o e v 

= 6 km.  With  these  profiles,  the  conductivities 
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7 

were  calculated,  taking  oe  = 10  rad/sec  as  the  electron 
gyro-f requency . From  the  calculated  conductivity  profiles, 
the  parameters  hc , c , h^,  and  were  determined  for 
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frequencies  of  50  Hertz  and  100  Hertz,  and  the  approximate 
eigenvalues  were  calculated  from  Eq.  (38) . A full-wave  cal- 
culation for  the  profiles  given  by  Eqs.  (52)  and  (53)  was 
carried  out  by  Dr.  William  Moler  of  Naval  Ocean  Systems  Center 
for  comparison  with  the  approximate  eigenvalues.  The  two  sets 
of  results  are  shown  in  Table  1.  The  phase  speeds  agree  to 
better  than  1%  and  the  attenuation  rates  to  within  0.2  decibels 
per  megameter,  which  is  quite  good. 

The  relative  phase  speeds  predicted  by  the  theory  of  Booker 
and  Lefeuvre  for  the  same  profile  are  approximately  0.95  at 
both  frequencies,  which  is  substantially  higher  than  the  full- 
wave  values.  We  have  not  attempted  to  calculate  their 
attenuation  rates,  which  requires  integration  over  the 
conductivity  profile  to  estimate  the  contribution  from  Joule 
heating  below  h^.  However,  an  important  assumption  on  which 
their  approximation  is  based  is  clearly  at  variance  with  our 
analysis.  They  calculate  the  fraction  of  energy  removed  from 
the  horizontal  flow  by  Joule  heating  at  a given  altitude,  and 
identify  the  attenuation  rate  with  the  average  value  of  this 
quantity  between  the  ground  and  h^ . This  identification 
relies  on  the  assumption  that  the  horizontal  flow  of  energy 
is  approximately  uniform  between  the  surface  of  the  earth  and 
the  level  of  reflection.  However,  the  lower  ionosphere  solu- 
tions presented  here  show  that,  while  the  horizontal  magnetic 
field  remains  essentially  constant  in  this  region,  the  ver- 
tical electric  field  falls  off  quite  rapidly  above  hQ.  To 
lowest  order  in  Co/hQ,  the  rate  of  horizontal  flow  of  energy  is 


Approximate-  and  Full-Wave  Propagation  Constants  for  the  Hypothetical 
Ionosphere  Described  in  Section  4. 


where  we  have  assumed  aQ  to  be  real.  For  an  exponential 
conductivity  profile,  the  horizontal  flow  rate  is  essentially 
constant  to  an  altitude  nQ,  above  which  it  falls  off  rapidly. 
The  Joule  heating  rate  is  given  by 


which  exhibits  a very  sharp  maximum  at  hQ,  and  becomes  very 
small  within  a scale  height  or  so  on  either  side  of  this 
altitude.  (These  points  are  illustrated  graphically  in  I.) 
For  the  profile  given  by  Eq.  (18)  , the  heating  and  flow  rates 
can  be  integrated  over  altitude  analytically.  Since  the 
integrands  fall  off  very  rapidly  above  h , little  error  is 
made  by  extending  the  upper  limit  to  infinity.  The  result  is 


/•°°  ? 

/ 0olEzl  dz 

*/o 

/«>  * 

Re(EzBx) dz 


= k \ Re(S)  , 
o 


which  agrees  exactly  with  the  horizontal  attenuation  rate  due 
to  heating  as  calculated  from  Eq.  (41) . 

Although  essentially  all  of  the  heating  dissipation  takes 
place  in  a narrow  altitude  region  around  hQ,  the  local 
attenuation  rate  (i.e.,  the  ratio  of  the  local  heating  rate 
to  the  horizontal  flow  rate  at  the  same  level)  is  an  increas- 
ing function  of  altitude.  Thus,  an  unweighted  average  of  the 
local  attenuation  rate  will  exceed  the  actual  horizontal 
attenuation  rate.  For  a given  profile,  the  difference 


23 
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I 

I I 

between  the  two  values  increases  with  decreasing  frequency 
due  to  the  lowering  of  the  altitude  hQ.  This  perhaps  accounts 
for  the  difficulty  experienced  by  Booker  and  Lefeuvre  (1978) 
in  reconciling  the  observability  of  the  nighttime  Schumann 
resonance  with  generally  accepted  nighttime  profiles. 


! 
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